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We introduce one dimensional sets to help describe and constrain the integral curves of an n 
dimensional dynamical system. These curves provide more information about the system than 
the zero-dimensional sets (fixed points) do. In fact, these curves pass through the fixed points. 
Connecting curves are introduced using two different but equivalent definitions, one from dynamical 
systems theory, the other from differential geometry. We describe how to compute these curves and 
illustrate their properties by showing the connecting curves for a number of dynamical systems. 
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I. INTRODUCTION 

Poincare proposed that the fixed points of a dynamical 
system could be used to provide some information about, 
or constraints on, the behavior of trajectories defined by 
a set of n nonlinear ordinary differential equations (a dy- 
namical system) PHI]. The fixed points of a dynamical 
system constitute its zero-dimensional invariant set. Un- 
fortunately, the fixed points provide only local informa- 
tion about the nature of the flow. 

Since that time many, including Andronov, Tikhonov, 
Levinson, Wasow, Cole, O'Malley and Fenichel, have fo- 
cused on higher dimensional invariant sets, in particular 
on ri — 1 dimensional invariant sets. In many instances 
these are slow invariant manifolds of singularly perturbed 
dynamical systems. These manifolds enable one to define 
the slow part of the evolution of the trajectory curve of 
such systems. Until now, it seems that, except for the 
works of [5] , no one has investigated the problem of one- 
dimensional sets which play a very important role in the 
structure of chaotic attractors by connecting their fixed 
points. The aim of this work is to define and present 
methods for constructing such one-dimensional sets. The 
sets that we construct are generally not trajectories that 
satisfy the equations of the dynamical system. 

The first attempt to study one-dimensional sets has 
been made in the context of Fluid Mechanics by Roth 
and Peikert [S]. The idea of transporting the concept 
of "vortex core curves" to the phase space of dynam- 
ical systems is due to one of us (R.G.), who applied 
it to three-dimensional dynamical systems and then to 
higher-dimensional dynamical systems. In the context 
of classical differential geometry, another of us (J.M.G.) 
called such curves connecting curves since they are one- 
dimensional sets that connect fixed points. 

In Sec. II we set terminology by introducing au- 
tonomous dynamical systems and define the velocity and 
acceleration vector fields in terms of the forcing equations 
for these systems. In Sec. Ill we introduce the idea of 
the vortex core curve through an eigenvalue- like equation 



derived from the condition that one of the eigendirections 
of the Jacobian of the velocity vector field is colinear with 
the velocity vector field. We show that this defines a one- 
dimensional curve in the phase space. In Sec. IV we in- 
troduce the idea of connecting curves from the viewpoint 
of differential geometry. These are defined by the locus of 
points where the curvature along a trajectory vanishes. 
In Sec. V we show that the two definitions are equivalent. 
In Sec. VI we describe three methods for computing the 
connecting curves for a dynamical system. Several appli- 
cations are described in Sec. VII including two models 
introduced by Rossler, two introduced by Lorenz, and 
a dynamical system with a high symmetry. The figures 
show clearly that the connecting curve plays an impor- 
tant role as an axis around which the flow rotates, which 
is why this curve is called the vortex core curve in hy- 
drodynamics. The figures also underline an observation 
made in [5] that this curve is at best an approximation to 
the curve around which the flow swirls. In the final Sec- 
tion we summarize our results and provide a pointer to 
visual representations of many other dynamical systems 
and their connecting curves. 



II. DYNAMICAL SYSTEMS 

We consider a system of differential equations de- 
fined in a compact E included in R™ with X = 

[n,i 2 ,..,i„]'e£cl": 



dX 
~dt 



(1) 



defines 



where Qt(X) = [fi(X), f 2 (X), f n (X)\ C 

a velocity vector field in E whose components fi are as- 
sumed to be continuous and infinitely differentiable with 
respect to all X{, i.e., are real- valued C°° functions (or C r 
for r sufficiently large) in E and which satisfy the assump- 
tions of the Cauchy-Lipschitz theorem [BJ. A solution of 
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this system is the parameterized trajectory curve or in- 
tegral curve X (t) whose values define the states of the 
dynamical system described by Eq. ([I]) . Since none of the 
components /j of the velocity vector field depends here 
explicitly on time, the system is said to be autonomous. 

As the vector function X (t) of the scalar variable t rep- 
resents the trajectory of a particle M, the total derivative 
of X (t) is the vector function (t) of the scalar variable 
t which represents the instantaneous velocity vector of M 
at the instant t, namely: 



dt 



(2) 



The instantaneous velocity vector ~$ (t) is tangent to the 
trajectory except at the fixed points, where it is zero. 

The time derivative of 1? (t) is the vector function j (t) 
that represents the instantaneous acceleration vector of 
M at the instant t 



7(<) 



dt 



(3) 



Since the functions fi are supposed to be sufficiently dif- 
ferentiable, the chain rule leads to the derivative in the 
sense of Frechet [7] : 



<93r dX 



dl_ 

dt dX dt 

as? 



(4) 



By noticing that ^fi is the functional Jacobian matrix 
J of the dynamical system ([!]), it follows from Eqs. ^ 
and Q that 



7 = J 



(5) 



This equation plays a very important role in the discus- 
sions below. 



III. DYNAMICAL SYSTEMS AND VORTEX 
CORE CURVES 

At a fixed point in phase space the eigenvectors of the 
Jacobian matrix define the local stable and unstable man- 
ifolds. At a general point in phase space the eigenvectors 
of the Jacobian with real eigenvalues define natural dis- 
placement directions. There may be points in the phase 
space where two eigenvalues form a complex conjugate 
pair and one is real, and the real eigenvector is parallel 
to the vector field that defines the flow. Under these con- 
ditions we expect that the flow in the neighborhood of 
such points swirls around the flow direction, much as air 
flow swirls around the core of a tornado. This parallel 
condition can be expressed in the coordinate-free form 



(6) 



The first equation is the mathematical statement of par- 
allelism; the second equation is a consequence of Eq.([5|. 

In coordinate form the eigenvalue condition can be 
written 



7i 



d_ , _ dfj dx 8 
dt dx q dt 



JW s = A/i l<i,s<3 (7) 



The condition that the acceleration field is proportional 
to the velocity field, ii = Ai; or /j = Xfi, can be repre- 
sented in the form 



A _ h _ [3 _ , 
A /2 h 



(8) 



The intersection of the surfaces defined by the first two 
equations defines a one-dimensional set in the phase 
space. This set is a smooth curve that passes through 
fixed points. Alternatively, the three equations define a 
one-dimensional set in the phase space augmented by the 
eigenvalue A: (xi, X2, £3, A). The projection of the one 
dimensional set from M 3+1 down to the phase space M 3 
defines the vortex core curve for the dynamical system. 

The arguments above are easily extended to define one- 
dimensional vortex core curves for n-dimensional dynam- 
ical systems. 

Eq. ^ has been used to try to identify the location of 
the "vortex core curve" [5] in hydrodynamic data. It is 
known that this equation provides a reasonable approxi- 
mation to the vortex core when nonlinearities are small 
but it becomes less useful as nonlinearities become more 
important [3]. 



IV. GEOMETRY AND CONNECTING CURVES 

The approach developed by Ginoux et al. [8HTU] uses 
Differential Geometry to study the metric properties of 
the trajectory curve, specifically, its curvature |1 IffTS] . A 
space curve is defined by a set of coordinates X (s), where 
s parameterizes the curve. Typically, s is taken as the 
arc length. When s is instead taken as a time parameter 
t, derivatives have a natural interpretation as velocity 
and acceleration vectors. The classical curvature along a 
trajectory is defined in terms of the velocity vector 1? (t) 
and acceleration vector 7 (t) by 



Si 



= Ki = 



7AV 



(9) 



Here 3? represents the radius of curvature. 
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We define connecting curves as the curves along which 
the curvature Ki is zero. 

Remark: Curvature measures the deviation of the 
curve from a straight line in the neighborhood of any 
of its points. The location of the points where the local 
curvature of the trajectory curve is null represents the 
location of the points of analytical inflection. 



A. First method 

In three dimensions, it is in principle possible to use 
two of the three equations 0y(X) = to express two 
of the three variables (x,y, z) in terms of the third, for 
example y = y(x), z — z(x). 



V. VORTEX CORE CURVES AND 
CONNECTING CURVE 

The dynamical condition Eq.([6]) that defines vortex 
core curves can be reexpressed as = 0. This is 

equivalent to the geometric condition Eq. ^ that defines 
connecting curves. As a result, the two definitions, one 
coming from dynamical systems theory, the other from 
differential geometry, are equivalent. 

Since the two definitions are equivalent, the conditions 
they provide for defining the connecting curve are also 
identical, as we now show. The vanishing conditions for 
the first curvature of the flow Eq. (|9| are 



= AT? JV?A T? 



<S> jA]f = 0^ ki = (10) 



By defining: 

012 = /1/2 
(c.f., Eq.(§) 



f>23 = /a/3 - fafs, 013 = /1/3 - /1/3 and 
/1/2 the third equality can be rewritten 
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I /2/3 — 72/3 — 
v = 6^l /1/3- A/ 3 = o 
A/2-A/2 = 




(ii) 



It can be proved that two of the three equations of this 
nonlinear system are equivalent and so this relation can 
be written as three subsystems: 




023 = 

012 =0 

013 = 

012 =0 

023 = 

013 = 



(12) 



By judiciously choosing one subsystem, say the first, 
we have another condition for defining the connecting 
curve, i.e. the intersection of two surfaces. 



^23 

f>12 



B. Second method 

If the dynamical system under consideration is of di- 
mension n the equation JV — XV represents a set 
of n equations in n + 1 variables: the n coordinates 
Xi,i — 1, ...,n and the eigenvalue A. These n equations 
define a one-dimensional set in the enlarged n + 1 di- 
mensional space. The projection of this one-dimensional 
curve into the n dimensional phase space is the connect- 
ing curve of the dynamical system. Since V — at the 
fixed points, all fixed points satisfy this equation and thus 
belong to the solution set. The method for construct- 
ing the parameterized version of the connecting curve 
involves writing down the n constraint equations for the 
71 + 1 variables (x, A), and eliminating all but one. 



C. Third method 

As previously observed, the problem for computing the 
connecting curve for three-dimensional dynamical sys- 
tems turns into the problem of computing the intersec- 
tion of two two-dimensional surfaces. A nice method 
for doing just this has been developed by Wilkinson 
P3]. We suppose that the intersection of two surfaces 
012 [x, y,z) = and 023 (x, y, z) = is parameterized by 
X (x (t) ,y (t) ,z (t)). The time derivative of the surface 

equation leads to V0ij (x, y, z) ■ X = 0. This means that 
X is perpendicular to both the gradients V0ij (x, y, z) 
which are the normal vectors to each surface. As long 
as these vectors are linearly independent for points on 

the intersection, then X is collinear to the cross product 
V0i2 (x, y, z) A V023 [x, y, z) 



X (t) = A (t) V0i2 (x, y, z) A V023 (x, y, z) 



(14) 



By rescaling the time t it is possible to set X(t) = 1. 
Then Eq. (14 1 simplifies to the form of an associated 



(13) dynamical system (A.D.S.): 



VI. CONNECTING CURVE COMPUTATION 

This kind of problem can not be solved analytically in 
the general case. As a result, three numerical approaches 
have been used to provide the connecting curve defined 
by the intersection of two surfaces, i.e., (13 1. 



dX(t) 
dt 



= V0i2 (x, y, z) A V023 (x, y, z) (15) 



These equations are generally different from, but related 
to, the original dynamical system equations. The curves 
defined by this equation are not heteroclinic trajectories 
of the original dynamical system. 
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Initial conditions for the (A.D.S.) are any point on the 
connecting curve, or any point belonging to the intersec- 
tion of both surfaces. Thus, the connecting curve may 
be defined as the trajectory, or integral, of the (A.D.S.). 
This method is useful as long as the gradients along the 
intersection remain nonzero and non-colinear. 



VII. APPLICATIONS 

In this Section we describe the connecting curves for 
three- and four-dimensional dynamical systems. 

A. Rossler model 



TABLE I: Coefficients of the fifth degree equation that de- 
fines the eigenvalue A in the expression for the curve along 
which the velocity and acceleration vectors are parallel for 
the Rossler dynamical system. 

D 5 = a 

D4 = 2a(c — a — x) 

D-i = ax 2 — 2acx + 4a 2 x — 4a 2 c + a 3 + c + 2a + ac 2 
Z>2 = — 2a 2 x 2 + x 2 - 2a 3 x — 2cx + 4a 2 cx — 4ax 

+ab + 2ac - 2a 2 c 2 + 2a 3 c + c 2 - 2a 2 
Di = a 3 x 2 + <ia 2 x-2a 3 cx-2a 2 b + a, + b + c-3a 2 c + a 3 c 2 
Do = x 2 — a 2 x 2 + 2a 2 cx — 2cx — 2ax + ac — a 2 c 2 + c 2 

-ab + a 3 b 



The flow equations for the Rossler attractor [TS] are 

/ i \ _> / h( x >y> z )\ ( -y-z \ 

V \y = 9f f2(x,y,z) = x + ay (16) 
\ z ) \f 3 (x,y,z) J \b + z{x-c)J 

where a, b and c are real parameters. The connecting 
curve for this dynamical system was computed using all 
three methods described in Sec. IV. The solution using 
the third method has been performed with Mathematica 
7 (files are available at: |http: / /ginoux.univ-tln.fr| . 

Of the three solution methods just described, the sec- 
ond leads to the simplest expressions for the connecting 
curve. The curve along which JV — XV depends on the 
three control parameters (a, b, c) and is parameterized by 
one of the three phase space coordinates. Choosing x as 
the phase space coordinate, the eigenvalue A satisfies a 
fifth degree equation 

5 

3=0 

The coefficients Dj are listed in Table I. At each fixed 
point, the value of A is the value of the real eigenvalue 
of the Jacobian matrix at that fixed point. The coordi- 
nates y and z are expressed as rational functions of x and 
A(a;; a, b, c). These rational expressions are 

— b — x + axle — x) + \x(x — c + a — A) 

y = 

a + (c — x){l — a 2 ) + Aa(c — x + A — a) 

(18) 

+b + x + (Ax + ab)(X - a) 
a + (c — x)(l — a 2 ) + Aa(c — x + A — a) 

The segment of the connecting curve between the fixed 
points (dots) is plotted for the Rossler attractor in Fig. 
[Tjfor control parameter values (a, 6, c) = (0.556, 2.0,4.0). 
Two projections are shown. Near the outer fixed point 
with repelling real eigendirection, this curve is a good 
approximation to a curve that defines the core of the 



tornado-like motion. However, as it moves toward the 
fixed point near the x-y plane and the nonlinearities in- 
crease in strength, it becomes a poorer and poorer ap- 
proximation of such a curve, even intersecting the at- 
tractor twice before joining the inner fixed point. This 
problem is apparent in the x-z projection. This result 
reinforces an observation made by Roth and Peikert that 
the "eigencurve" , i.e. the connecting curve defined by 
JV = XV, is a good approximation to the vortex core 
curve in regions where the nonlinearities are weak, but 
not where the nonlinearities become strong [S]. 



B. Lorenz model 

The purpose of the model established by Edward 
Lorenz [16] was initially to analyze the unpredictable 
behavior of weather. After having developed non-linear 
partial differential equations starting from the thermal 
equation and Navier-Stokes equations, Lorenz truncated 
them to retain only three modes. The most widespread 
form of the Lorenz model is as follows: 




*\ / A (x,y,z) \ I a(y-x) \ 

y = Of h(x,y,z) J = I Rx — y — xz\ (19) 

z ) \h{x,y,z) ) \ -bz + xy J 



where a, R and b are real parameters. Once again, the 
connecting curves were computed using all three meth- 
ods described in Sec. IV. The calculation using the 
third method was performed with Mathematica 7 (Files 
are available at: http://ginoux.univ-tln.fr). All methods 
gave the same curves. 

Three connecting curves pass through the saddle at the 
origin: one corresponding to each of the three eigendirec- 
tions with real eigenvalues. The simplest of these curves 
is the z-axis, which is simple to compute by hand. This 
particular curve is a trajectory of the Lorenz model. A 
second heads off to z — > —00 and has little effect on the 
attractor. The third connecting curve passes through all 
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FIG. 1: Connecting curve of the Rossler model. The curve 
intersects the attractor twice, as seen in the x-z projection. 
Parameter values: (a,b,c) — (0.556,2,4). 



three fixed points. This curve is shown in Fig. [2] in both 
the x-y and y-z projections for (R,a,b) = (28,10,8/3). 
When R is increased, the return flow from one side of the 
attractor to the other exhibits a fold and the connecting 
curve intersects the attractor at the fold. This reflects a 
similar property shown by the Rossler equations. 

The connecting curves present additional constraints 
on the structure of the Lorenz attractor above and be- 
yond those implied by the location and stability of the 
fixed points. Specifically, the flow spirals around and 
away from the connecting curve that passes through the 
two foci. In addition, the z axis also provides some struc- 
ture on this flow, as the flow also always passes in the 
same direction around this axis |17j . 



C. Lorenz model of 1984 

In 1984 Lorenz proposed a global atmospheric circula- 
tion model in truncated form [IB] . The model consists of 
three ordinary differential equations: 



FIG. 2: Connecting curve of the Lorenz model. One nontrivial 
connecting curve heads off to z — > — oo and has little effect on 
the structure of the attractor. The other nontrivial connecting 
curve connects all three fixed points, and is plotted extending 
through the foci. The third connecting curve is the z axis. 
Parameter values: (R, a, b) = (28, 10, 8/3). 



fi (x,y,z) 
h {x,y,z) 
h {x,y,z) 

' —y 2 - z 2 - a(x 
-y + xy — bxz 
bxy + xz — z 




(20) 



In this model the variable x represents the strength of the 
globally circling westerly wind current and also the tem- 
perature gradient towards the pole. Heat is transported 
poleward by a chain of large scale eddies. The strength 
of this heat transport is represented by the two variables 
x and y, which are in quadrature. The control parame- 
ters aF and G represent thermal forcing. The parameter 
b describes the strength of displacement of the eddies by 
the westerly current. 

In Fig. [3] we show two projections of this attractor for 
control parameters (a,b,F,G) = (1/4,4,8,1) as well as 
the connecting curve. For this set of parameter values 
there are three fixed points, only one of which is real at 



G 



(x,y,z) = (7.996,-0.00653,0.0298). It is clear that the 
connecting curve goes through the hole in the middle of 
the attractor, and that the attractor winds around part 
of the connecting curve where most of the bending and 
folding of the attractor occurs. The connecting curve in 
the x-y projection passes through the fixed point off scale 
to the right. 








N 




FIG. 3: Strange attractor generated by the Lorenz global cir- 
culation model of 1984. The connecting curve threads through 
the inside of the attractor, and is caressed by the attractor 
where the stretching and folding is most pronounced. Param- 
eter values: (a,b,F,G) = (1/4,4,8, 1). 



fx 



z 
w 



/ h (x, y, z, w) 

f 2 (x, y, z, w) 
h (x, y, z, w) 




(21) 



Here the state variables are (x, y, z, w) and the control pa- 
rameters are (a, b, c, d). The connecting curve was com- 
puted using methods 1 and 2 of Sec. IV. The first method 
gave very complicated results. Method 2 gave simpler 
results when the coordinate z was used to express the 
behavior of the remaining four variables. The eigenvalue 
A was expressed as the root of a seventh degree poly- 
nomial equation whose coefficients were functions of the 
four control parameters (a,b,c,d) and z. The remaining 
three coordinates were rational functions of small degree 
in the variables z and \(z; a, b, c, d). Two projections 
of the hyperchaotic attractor and the connecting curve 
are shown in Fig. [4j The computation was carried out 
for (a,b,c,d) = (1/4,3,1/2,1/20). The fixed points are 
shown as large dots along the connecting curve. It is 
clear from this figure that the connecting curve provides 
information about the structure of the attractor, as the 
flow in the attractor swirls around the connecting curve. 



E. Thomas Model 



Thomas proposed the following model of a feedback 
circuit with a high degree of symmetry [20] : 



h (x,y,z) 
h (x, y, z) 
h (x, y, z) 



ay 

az - 
ax 



(22) 



D. Rossler model of hyperchaos 

Rossler proposed a simple four-dimensional model in 
1979 to study hyperchaotic behavior [TJ5]. This model is 



This set of equations exhibits the six-fold rotation- 
reflection symmetry Sq about the (1,1,1) axis. The 
symmetry generator is a rotation about this axis by 
27r/6 radians followed by a reflection in the plane 
perpendicular to the axis. The origin is always a 
fixed point and, for a — b > 0, there are two on- 
axis fixed points at x = y = z = ±\/b — a. For 
(a, b) = (1.1,0.3) there are 24 additional off-axis 
fixed points. These fall into four sets of symmetry- 
related fixed points (sextuplets). One point in each 
sextuplet is (0.085,1.037,0.309), (0.250,1.013,0.865), 
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FIG. 4: Hyperchaotic attractor generated by the 1979 Rossler 
model for hyperchaos. Parameter values: (a, b, c, d) = 
(1/4,3,1/2,1/20). 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



y 

FIG. 5: Connecting curves for the Thomas attractor. One 
connecting curve is the symmetry axis x — y — z. The re- 
maining connecting curves exhibit the six-fold symmetry of 
the system and seek out the holes in the attractor. Parame- 
ter values: (a,b) = (1.1,0.3). 



(0.364,-1.095,1.175), (1.146,-1.180,-0.816). The re- 
maining points in a multiplet are obtained by cyclic per- 
mutation of these coordinates: (u,v,w) — > (w,u,v) — > 
(y,w,u) and inversion in the origin (u, v, w) — > 
(—u,—v,—w). The chaotic attractor for this dynamical 
system is shown in Fig. [5j along with the symmetry- 
related connecting curves and the 27 fixed points. One 
of the connecting curves is the rotation axis. This is an 
invariant set that connects the three on-axis fixed points. 
It therefore cannot intersect the attractor. In fact, this 
set has the same properties as the z-axis does for the 
Lorenz attractor of 1963 [17]. The remaining connect- 
ing curves trace out the holes in the attractor. In this 
sense they provide additional constraints on the struc- 
ture of the attractor over and above those provided by 
the spectrum of fixed points. 



VIII. DISCUSSION 

In this work we go beyond the zero-dimensional in- 
variant sets (fixed points) that serve to a limited extent 
to define the structure of an attracting set of a dynam- 
ical system. We have introduced a curve that we call 



a connecting curve, since it passes through fixed points 
of an autonomous dynamical system. We have defined 
this curve in two different ways: dynamically and kinc- 
matically. It is defined a vortex core curve dynamically 
through an eigenvalue-like equation JV = XV, where 
V(x) is the velocity vector field defining the dynamical 
system and Jy = dVi/dxj is its Jacobian. We have de- 
fined a connecting curve kinematically as the locus of 
points in the phase space where the principal curvature 
is zero. These two definitions are equivalent. 

Three methods were introduced for constructing this 
curve for autonomous dynamical systems. They were 
applied to the standard Rossler and Lorenz attractors, 
where their behavior with respect to the attractors is 
shown in Figs. 1 and 2. In the figures shown for these 
attractors, it is clear that the flow rotates around the con- 
necting curves, which therefore help to define the struc- 
ture of the attractor. The connecting curves were also 
constructed for a later Lorenz model, the global atmo- 
spheric circulation model of 1984, and for a later model 
introduced by Rossler to study chaotic behavior in four 
dimensional phase spaces. Finally, a multiplicity of con- 
necting curves was computed for an attractor with a 



8 



high degree of symmetry, the Thomas attractor. This 
is shown in Fig. 5. The flows shown in Figs. 3, 4, and 
5 are clearly organized by their connecting curves. In 
this sense the connecting curve provide additional im- 
portant information about the structure of an attractor, 
over and above that provided by the number, nature, and 
distribution of the fixed points. A number of other con- 
necting curves have been computed, and can be seen at 
http: / /www. physics, drexel.edu/~tim/programs/ 
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